IMSc-96/09/28 



On a Static solution of Einstein Equations with 
incoming and outgoing radiation 



G. DateQ 



Abstract 

Einstein equations with T^ u = k^k u + £^£ u where k, I are null are 
considered with spherical symmetry and staticity. The solution has 
naked singularity and is not asymptotically flat. However, it may be 
interpreted as an envelope for any static spherical body making it 
more massive. Such an interpretation and some of its implications are 
detailed. 
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Introduction 



Any astrophysical body has a temperature and thus is a source of out- 
going radiation. Equally well, every such body also receives an incoming 
radiation eg. the cosmic background radiation. Depending upon the respec- 
tive temperatures there will either be a net outgoing or incoming flux of 
radiation. It is then conceivable that the two rates exactly match and one 
reaches an equilibrium situation. In such a situation, the stress tensor in the 
vicinity of such a body may be taken to be of the form 

T^v = kfjky + £„£ v , k 2 = f = , k.l > (1) 

When the rates are unequal either the k or the £ term may be taken to 
be dominating and one essentially gets the Vaidya solution (non stationary) 
[[!]], or collapsing null fluid shell case. However when the rates are precisely 
matched, both terms are important and one can look for a static solution. 
As a first step in this direction of course one can consider the simpler case 
of spherically symmetric solution. 

In the present work we consider such a solution. Some of the salient 
features of the solution are the following. 

Let r denote the usual Schwarzschild radial coordinate. 

1. As r — > the solution has a curvature singularity which is naked (i.e. 
no event horizon). 

2. As r — > oo, the metric components go as ln(r) and thus the solution is 
not asymptotically flat. 

3. These two features make it difficult to interpret the solution physi- 
cally. However, one can consider the solution to be valid for R < 
r < R range. At R one can match the solution for a typical interior 
Schwarzschild solution while at R one can match it with the standard 
exterior Schwarzschild space-time. The asymptotics of the solution 
permit such a (C°) matching. If M denotes the mass of the interior 
solution and M denote the mass indicated by the matching at R, thus 
M > M and therefore the mass measured (deduced) from r >> R is 
larger than M. One can consider the matching at R just outside a 
black hole or even with a negative mass Schwarzschild solution, but M 
is always positive for sufficiently large R. 
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The paper is organised as follows. 

Section 1 contains basic equations which are straight forward to derive 
In Section 2, we present analysis of qualitative features of the solution. The 
asymptotics are discussed and numerical solutions are presented corroborat- 
ing the qualitative analysis. In Section 3 we discuss some of the possible 
matchings and summarise our conclusions. The appendix contains a few de- 
tails of general spherically symmetric and static non-empty space time. A 
some what similar exact solution is also presented and the details of match- 
ings are discussed. 

Section 1: Basic equations 

The basic equations are (Signature H ) 

Rf,u - ^Rg^u = 8nT^ (2) 

Tpv = PKK + <tWu , k 2 = I 2 = , kl > . (3) 

where vector fields represent massless radiation outgoing and in- 

coming respectively. 

Clearly g^T^ = and therefore R term can be dropped. Using spherical 
symmetry and staticity we write 

ds 2 = F(r)dt 2 -G(r)dr 2 -r 2 (d9 2 + sin 2 9d<p 2 ) (4) 
The killing vectors are : 



t| — > generating time translations ; 
7) = + where, 

(5) 

£ e = a sin + (3 cos 

^ = 7 + cot 9 (a cos t(p — (3 sin 0), 



generating isometries implied by spherical symmetry. 
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It follows immediately that 

C ( T^ = <-> £^ = = (6) 

Combining with k 2 = I 2 = one gets 

= fc°(l,Jg,0,0) , €" = *°(1,-Vg,0,0) 



A; , £° > , 2k°£°F > W 

Here /x = 0,1,2,3 <-> t,r,Q,<j) respectively. Since are arbitrary at 

this stage, we absorb p and a in fc° , £° respectively. 

We can actually solve for k° and £° using the conservation equations. 

V„ T^ u = =^ 

(jfe"V ■ k + k ■ VJfe") + (rV-£ + £-V^) = 0. ^ 

Evaluating the covariant derivatives etc. shows that each bracket is sepa- 
rately zero implying that k^ and £^ integral curves are geodesies (non affinely 
parametrised). Furthermore, one gets 

B B 

k° = — ^ , £° = — -, where B± are constants. (9) 
rF rF 

Therefore, 

** = % t 1 '/!' 0,0 )'" = § i i --v5' o '°i <io) 

The nonvanishing components of then are 

T - B + + B - T -°T T - F ( B -- B2 A (ll) 

For spherically symmetric, static metric i?oi is zero an d therefore = 
B 2 _ and since k°,£° are both positive we take B + = B_ = £> . Then, 

k° = £° = ^- , A = 2£? 2 (13) 



Note that in the geometrised units we are using T has dimensions of 
(length) _2 and thus A is dimensionless. Since T 00 ~^we see both that we 
cannot have asymptotic flatness and that there will be a curvature singularity 
at r = 0. 

It is straightforward to compute R^ u . Out of the 10 Einstein eqns. the 
eqns for (fj,u) = 02, 03, 12, 13, 23 are identically satisfied. The 01 equation has 
already been used to set B + = B_. The 22 and 33 equations are identical. 
This leaves us with 3 nontrivial equations. Setting / = \f(F) g = \f{G) , 
the equations are: 



(00) 

(11) 
(33) 



-r 1 K) + 



-r 



-if 

9 2 



+ 



2r- x 4 
g 2 

2r~ 1 4 

9 







(14) 



Remarks 

1. A = gives the standard Schwarzschild case. The 33 eqn is independent 
of A and is invariant under constant rescaling of /. 

2. Under / — > A/ the 00 and 11 equations retain their form but with 
A -> A/\ 2 

3. All the equations are also invariant under r — > Ar and therefore there 
is no intrinsic scale available at the level of the equations. This has to 
be provided by physical boundary conditions. This is of course true 
for Schwarzschild case as well (indeed whenever the stress tensor is 
traceless and has no dimensionful parameters). 

4. By suitable combinations at 00 and 11 equations we get one equation 
which contains A dependent term while another one which does not 
contain explicit A dependence. The A independent combination is a 
second order differential equation while A dependent one is a first order 
equation. The 33 equation is 1st order. 

It is straightforward to verify that the second order equation is automat- 
ically satisfied if the two first order equations are satisfied. 
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Defining r = /io^, Ho an arbitrary scale and A = 8nA, we write the 
equations in terms of F and G, as : 

dF 

^ = XG + F(G - 1) (15) 

A^-G(G-l) (16) 

2d' 2 F fid V _ dF dG 2dG 
F~d~e ~ \~Fdi J ~F^Gdt, ~ G~dJ ~ ( ' 

The last equation is the second order equation which is identically satisfied 
if the first two equations hold. Thus the basic equations to be solved are the 
first two equations. 

To summarise : 




ds 2 = Fdt 2 - Gdr 2 - rW 

fc" = fe°(l, > /g,0,0) 
= ^(1,-^/1,0,0) 

r = /x e« , / <-> ^ 

F = AG + F(G - 1) 
G' = A^-G(G-l) 

The case A = gives the Schwarzschild solution. 

For A 7^ 0, we set F = A$. The equations for G and <£> then have no A 
dependence and are: 

$' = G + $(G-1) (A) 

G' = ^-G(G-1) (5) 
In the next section we will analyse these equations. 
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Section 2: Analysis of the equations 



$ > , G > 

$' = g + $(G-1) , G" = ^-G(G-l) 1 J 

Therefore, 

n $ + 

Eliminating G from the second equation, gives a second order equation for 
<3>, namely, 

$"{$($ + 1)} + $'{$($- 2)} -$ /2 -2$ 2 = 0. (19) 
This can be integrated once to give, 

$'(1 + -) + $- 21n($) -2£ = C (20) 

or, 

$' = (C-$ + 21n($) + 20( ¥ ^ T ) (21) 
Substituting in the expression for G, we get the exact equations: 



G = [C + 1 + 2 ln($) + 2f ] (solves G in term of $) 



(22) 



$' = (^Jl) [C - $ + 2 ln($) + 2f] (1st order equation for $) 



The first integral has given us one constant of integrations. 
Remarks: 

1. G = 0, $' = — $ is an exact solution. This follows from both the 
original equations for <3> and G and from the above equations. However, this 
is not an acceptable solution. 

2. If $' = then 

C + 2 ln($) + 2£ = $ (23) 



$ + 1/ L J V$ + l 



-$' + -$' + 2 



(24) 



Therefore, 

2$ 

$'V=o = + - — - >0 (25) 
$ + 1 
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Thus $ has at the most one minimum. The minimum is determined by 



$(0-21n($(0) =C + 2i (26) 
Therefore if <&(£) = a then a — 2 ln(a) = C + 2£ determines a given C and 

i 

By adjusting yUo, we can always choose £ to be zero or any value. 
At the minimum of $, $ = C + 2 ln(<l>) + 2£ which implies: 

G = (*TT? I$+11 = ¥TT < 1 

If $ — > at any finite £ then $' — > as £ — > £ and G — > as well. Since 
at finite £ we have everything regular, the det g ^ and therefore $ and G 
both can not be allowed to go to zero. 

Thus, there can be no event horizon at any finite £ , and $ and G are 
necessarily > V finite £. 

Now consider the asymptotics: 



Observe that $ can not be oscillatory since it can have at the most one 
extremum. Hence it is either bounded or unbounded as £ — > ±oo. 

If $ has a finite, non-zero limit, then the equation implies that $' — > ±oo 
which is absurd. Therefore $ either vanishes or diverges to oo. 

If $ — > , then we can approximate the equation for $' as: 

$' ~2 $ (£ + ln($)) =>■ 

ln($) ~ #e ^-I-£. 

This will go to — oo provided # is zero and £ — > +oo. Clearly then as 
£ — > — oo, $ must diverge to infinity. 

As £ — > +oo, $' — > which in turn implies that G — > as well. But 
= 2G 2 and therefore $G increases monotonically and hence can noi 
vanish. Thus $ it must diverge as £ — > +oo as well. 

To summarise, $ — > +oo as £ — > ±oo must hold which in turn implies 
that <£> must have a minimum. 

Consider approximate solution as $ — > oo(£ — > ±oo) . 

Let 

77 = $ - 2 ln($) - 2£ - C. (27) 
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Therefore 



rf = |^f(-^)- 2 where$ = Hv) (28) 
Expanding in powers of 1/$, 

= -^-Iw-^ + ^-Ti-)- 2 v *>!- (») 

For $ >> 1, 77' ss -77 - 2 =>- 

= £> e -*-2 =>■ 

$-21n($) = C-2 + 2£ + J De _$ (30) 

and this is consistent with $ >> 1 provided either £ — > +00 or £ — > —00. 

The corresponding asymptotic behaviour for G can be deduced from the 
behaviour of $. The leading behaviours are given below. 

£ — > cx) : $ k 2( 

: G « 1 + i 

£ — > -00 : $ ft De~^ 

: G ft £ (G + l + 21n(D)) 

Remark: The £ — > — 00 behaviour shows that $G — > constant. We may 
choose this constant to be 1 by choosing £> = e~ c ^ 2 . The asymptotic form 
then resembles the form for a negative mass Schwarzschild solution. This is 
not surprising since as £ ^ — 00, G — > and $ — > 00 , the first terms in 
the basic equations become negligible and the equations approximate to the 
standard Schwarzschild equations. 

To summarise : The equations can be integrated once exactly to give 

1. 

$' = (^)[C< + 2£-$ + 21n($)] 
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2. Regularity at finite £ =>- $, G > V finite £. Therefore no event horizon 
is possible. 

3. $ can have at the most 1 extremum which must be a minimum. 

4. $ — > oo as £ — > ±oo and therefore $ does have the minimum. 

5. G" > as ^ — > — oo while G' < as £ — > oo and therefore G" must 
vanish for some finite £. G then /ias a unique maximum. 

This qualitative picture is corroborated by numerical integration of the 
$ — G equations as shown in the figure below. The specific initial values are 
chosen for convenience only. 




Figure 1: The solution has G(0) = 1.25, G'(0) = 0.0 



Section 3: Interpretation of the solution 

We have four freely specifyable parameters: the arbitrary scale /i , the 
constant of integration A coming from the conservation equation, the constant 
C and a = $(0) (say). After the substitution F = A$, the A drops out of all 
differential equations. It still appears in the time components of the vector 
field k, I and the 00 component of the stress tensor. It also appears in the 
metric as a coefficient of the dt 2 term. By rescaling time we can remove it 
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from the metric and from the vector fields k, £. Clearly then its value cannot 
have any physical meaning. In effect we take A = 1. 

The scale p Q on the other hand appears non trivially in the following 
sense. If we take any stationary observer with four velocity w M = £ M / y/W, 
then the energy density measured in his/her rest frame is given by, 

p u (r) = T^u v 

~ 87rr 2 $(r) (31) 
J_ 1 e~ 2 £ 
Stt ^ *(£) 

How could a scale be chosen? As noted earlier, our solution has a curva- 
ture singularity as r — > 0. So a natural approach is to consider the solution to 
be valid for r > R for some R. This R then provides a natural scale (po = R). 
We also noted earlier that the £ going to — oo behaviour resembles that of 
a negative mass (— m) Schwarzschild solution. The constant D then equals 
2m/ Ho. The parameter 2m then provides a natural scale. In either of the 
cases, a natural scale p,o can be chosen. It remains now to choose C and a. 

If matching with negative mass Schwarzschild solution is considered then 
the choice po = 2m gives D = 1 or C = 0. The constant a is left uncon- 
strained (> 0). 

A more "realistic" matching is to choose an R, the radius of some physical 
body and match our solution with an interior Schwarzschild solution (perfect 
fluid case for instance) . For an interior solution, the function G is expressed 
in terms of a mass function M(r) 

M(r) = An f p(r')(r'fdr' 
Jo 

and the equations are integrated (usually numerically). The $ function is 
trivially determined once p(r) is determined (The equation of state gives the 
pressure P(r) ). The matching is minimally required to have F, G and F' to 
be continuous across the matching surfaces (See the appendix for details). 
The continuity of $ across R can always be ensured trivially since for the 
interior F there is a constant of integration which can always be adjusted. 

A physical body provides the following physical data namely, the physical 
radius R which gives po] the physical mass M (in the absence of the radia- 
tion shell) which gives G(R) = (1 — 2M/R) ; and the radiation density 
measured by a stationary observer in his/her rest frame p u {R) which gives 
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As discussed in the appendix, p u {R) = P(R) because of continuity of 
$' . Thus all the data needed for specifying a particular solution is available. 
Thus for interior matching, we take: 

b = G{0) = 



Ml 
. R 



» - •<«) = <32) 
Given a and 6, the constant C is given by, 

C = &(Q + 1)2 -l-21n(q) (33) 

This gives a method of choosing the constants of integration in a given 
physical context, thus determining the solution appropriate for the context. 

The solution so determined is to be evolved up to some finite £ as we 
do not have asymptotic flatness. At this point we would like to match our 
solution to an exterior Schwarzschild solution. As discussed in the appendix, 
it is not possible to do so while maintaining the continuity of <£>' and a "reg- 
ularising" thin layer must be added on. This can be done. Suppressing the 
"thicknes" of the regularising layer we see that the continuity of G provides 
us with an M. However, since $(£) is not equal to G~ 1 (£) the exterior F 
function will go to $(£) as £ goes to infinity. The mass given by the Komer 
integral will then have a normalization such that this mass is given by the 
M. 

Thus for exterior matching at £, we set 

1 R e 



It follows then, 



M 



M(R) 



g(0 - 1 G(0) 
G(0)-1 G(0 



e* (35) 



Mass of such a body will be larger by factors, from the mass it would 
have had in the absence at the radiation shell, (equivalently from the mass 
determined from the interior dynamics). 

To get a feel, let us put in some numbers. Let At denote the energy 
density of the background radiation at temperature T. It is given by, 



A T ~ KT 15 x T 4 ergs/cm 3 
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The p u {0 on the other hand is given in conventional units by, 



c 4 ,,1 w e~ 2 « 



^Newton oTTU <±> 



Therefore, 
Or, 

Putting £ = gives, 



r 4 



1 , 



e* ~ (10 31 )( 



10-) ( 



i , 



' i? 2 T 4 

For white dwarf (say), b is typically about 1 + 10~ 4 0. Typical white dwarf 
radius is about 10 9 cm. An astrophysical body such as a white dwarf could 
not be expected to have been formed in the earlier epoches and thus the 
background temperature can not be larger that about 10 4 . a is then about 
10 28 ! 

Notice that in this case a(b — 1) ^> b. As £ is increased a increases and b 
decreases, relatively slowly, maintaining the inequality. Thus the first terms 
in both of the basic equations are negligible. But then the equations ap- 
proximate the usual Schwarzschild case and M(£) read off from G will be 
essentially a constant i.e. M/M will be very close to 1. Numerical corrobo- 
ration of this is shown in the figure below. 

It is also clear that to get significant deviations from the usual Schwarzschild 
case one must have a(b — l)/b to be comparable to 1 or less than 1. In any as- 
trophysical context (excluding black holes), the T would be about the same 
order while b — 1 continues to be not too small. Only way then to get a 
deviation is to reduce a i.e. increase R. 

Indeed if we take a spherical galaxy to be the inner body then R ~ 10 22 
cms. (10 5 light years), the mass M is about 10 12 solar mass giving b— 1 ~ 10 -4 
- 1CT 5 . a is about 100 and significant deviations from the usual Schwarzschild 
solution can be expected. The figure below corroborates this expectation. 
The ratio M/M is several orders of magnitudes bigger than 1 even for R/R 
slightly larger than 1 ! The radiation shell will of course be dark since it 
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5 



1 



Figure 2: The solution has $(0) = 10 10 ,G(0) = 1 + 1(T 6 

will be merged with the background radiation. The shell then seems to be a 
candidate for dark matter at least in some cases. 

What about the core being a black hole? Clearly we can not match our 
solution at the horizon but we could try matching outside the horizon. A 
black hole can provide an out going flux only via the Hawking mechanism. 
Since we are using stationary observers to get the value of a, the energy den- 
sity should be that corresponding to the local temperature. The temperature 
is then given by, 



For R not too large {(5 of the order of 10, say) deviations are possible only 
for extremely light black holes! 

Remark Quite apart from these numbers, the precise matching of the 
two rates is fragile though. As the universe expands, the incoming rate will 
decrease and the black hole will begin loosing its mass, thereby increasing its 
out going rate and preventing return to equilibrium. The net result of the 
earlier equilibrium is perhaps to delay the evapouration process. 



Setting a = M/M Q and (3 



a(b 
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1 1.02 1.04 1.06 1.08 1.1 

Figure 3: The solution has $(0) = 10 2 , G(0) = 1 + 1(T 5 

The similarity of our solution to a negative mass Schwarzschild solution 
for r close to zero suggests a purely speculative possibility of taking the core 
to be a "negative mass" body. Of course no such body is known! If at all it 
"exists" one could only imagine a quantum origin. All the scales may then 
be taken to be Planckian. a then is of the order of 1. Numerical solution 
then indicates that the M becomes positive for R/R greater than about 2 - 
10. 

Though spherical symmetry, staticity and the particular form of the stress 
tensor are obvious idealizations, one can still observe the following: 

1. The stress tensor satisfies all the usual energy conditions and as such is 
a physically possible/admissible one. The Einstein equations then lead 
us to a solution discussed above. 

2. Apart from the role of providing incoming radiation, the background 
radiation or its temperature appears explicitly in providing one of the 
constants of integration, $(0). The solution then seems best interpreted 
as a radiation shell near the surface of a spherical static body. This shell 
however contribute to the mass significantly only for bodies with sizes 
on the galactic scale. The effective darkness of such bodies indicates a 
possibility for dark matter at least in some cases. 
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Appendix 



We collect here a few details of the most general equations for static, 
spherically symmetric non-empty space times which are useful for discussion 
of the matchings involved. 

Define the orthonormal set of vectors: 



^(1,0,0,0) ; ef = ^(0,1,0,0); 
e£ = ±(0,0,1,0) ; e£ = -1^(0,0,0,1). [ } 

Any T^ v can then be expressed as, 

= p ab e%el, p ab = p ba . (2) 

Spherical symmetry and staticity (i? i — 0) implies that, 

p ah = diag(p (r),p 1 (r),p 2 (r),p 3 (r)), (3) 

with p 3 (r) = p 2 (r). 

The conservation equations give a single equation: 

d Pi , 2 (Pi -P2) (p + pi) dF = 

dr r 2F dr 1 ' 



Some special cases are: 

Perfect fluid 
Reissner- Nordstrom 
Present case 



Po = p , p! = p 2 = P 

Po = ?r , Po = P2 = -Pi (5) 
p 2 = , pi = p = p 
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The Einstein equations can be organised as explained in section 1 to get 
a set of three first order differential equations as: 



dF 

dr 



F(G - 1) + (87ir 2 Fpi)G .... (a) 
f = -G(G-l) + (87rr 2 Fp )f .... (6) (6) 
£ = 2( P2 - Pl )-^f •••• (c) 



Defining, 



(7 a = 87rr 2 Fp a a = 0,1, 2 
and using the dimensionless variable £ = ln(r//io), we get: 

: F(G — 1) + o\G .... (a) 



dF 
d( 



dG _n(n _ i ^ i ~ G 



G(G-1) + <t <£ .... (6) (7) 

d i£ = 2^2 -^f .... (c) 

Mathematically we have an underdetermined system of equations with <t 
and o"2 (say) as freely specifiable function. Physically of course the a a 's are 
to be determined by the dynamics of the matter constituents eg. Maxwell 
equations for the Reissner-Nordstrom case, equation of state for the perfect 
fluid case and modeling in our case. 

Usually one notes that the (6b) equation does not involve F and solves 
this equation in terms of, 

M(r) = 4vr | p Q (r')(r') 2 dr' , G = (1 - ^^V 1 

The (6c) equation which involves only p a 's is then solved and (6a) equa- 
tion can then be trivially integrated. 

Since in our case equation (6c or 7c) are trivially integrated we were able 
to reduce the remaining equations to a single first order differential equation. 

In fact following our way of organising the equations one can construct 
the following exact solution . 

Taking o\ = o"o = a and 02 = —Ae~^ gives us a = B + 2Ae~^ and allows 
the remaing two equations to be reduced to a single first order equation for 
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F. A = reproduces our equation (22) of section 2. However for B = 
and A > (to satisfy energy conditions) further integration gives an exact 
solution (F = A$, A > 0) : 



$ + 2e _ * ln($) 
G 



= C + (£-20e-« 



($+2e-«) 2 
e -3e 



c 



C > 



(8) 



P/2 



This solution again has naked singularity and although p falls off faster 
than before, it is still not fast enough to get asymptotic flatness. 

Now we address the issue of matching our solution on the interior to a 
physical body and on the exterior to an exterior Schwarzschild solution. 

It is useful to note that the energy flux, as implicit in the Komer integral, 
depends on derivative of F. The extrinsic curvature for the hypersurfaces 
t — i constant and r = constant also depend on the derivative of F apart from 
on F, G, r etc. None depend on the derivative of G though. If the matching 
are to ensure continuity of any (or all) of these, then one must demand 
continuity of F, G and F' across the matching spheres. The equations 6a 
(or 7a) show immediately that such a matching can not be done on the 
exterior because p\ is not continuous across R . On the interior though such 
a matching is possible. Note that the matching requires continuity of p\ only 
and NOT of p . 

For an interior body described by perfect fluid stress tensor we have, 



Putting a = F(0), b = G(0) d = F'(0) and noting that b = (1 - 
2M(0)/R)- 1 we see that 




dF 



2G(0) 



M(0) + 47ri? 3 P(0) 
R 



d = a(b - 1) + 8irR 2 abP(0) ( for interior body) 



while for our solution 



d = a(b - 1) + b 



Thus the matching implies, 



1 



a = 



8nR 2 P(0) 
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By comparing with equation (31) of section 3, we see that P(0) = p«(0). 
This fixes the constant of integration for the interior solution as well as 
provides an initial condition for our equations. 

In the absence of the radiation shell, R is determined precisely by de- 
manding P(R) = 0. We can not do so. So our choice must be some what 
smaller than the usual value for R. 

Recall from section 3 that we chose a by equating the energy density 
measured by a stationary observer to the background energy density. In 
the conventional units, p«(0) (and hence P(0)) is about 10~ 15 x T 4 which is 
small enough so that the value of R will be very close to that determined in 
the absence of the radiation shell. So for estimate purposes we can use our 
earlier estimates. The solution thus determined is essentially the same as in 
the section 3. 

The exterior matching is to be examined now because the estimate of M 
depends crucially on this. 

We can modify the matching at £ used in section 3 by adding a thin 
"regularising layer" which will match with our solution at R and match with 
exterior Schwarzschild solution at slightly farther away. This matching of 
course is to have continuity of F' as well. 

To describe the thin layer, we consider the general equation (7). We 
retain cr 2 = condition but allow o\ ^ a . We will choose o\ suitably and 
solve for <7o using eqn(7c), i.e. 

0-0 = 0-1- 2o{ — . 

We choose a\ such that: 

<ri(0 = 1 ; <7i(£+2e) = 0; 

= ; a[(t + 2e) = 0. W 

This implies that 0"o(O = 1 and <t (£ + 2e) = 0. By choosing U\ to be 
monotonically decreasing we can ensure tha energy conditions are satisfied. 
A simple choice is: 

= le-frg-a'K *"^" 6 ) 3 - 3( e ~^~ e ) + 2} (10) 

By taking v large one can control the decrease in o\ while by taking e 
small one can make the layer thin. At £ + 2e one can match with the exterior 
Schwarzschild solution and read off mass from the value of G. 
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Numerical exploration of this thin layer shows that the qualitative con- 
clusions derived in the section 3 do not change. In particular, for stellar 
scales the radiation shell contributes negligibly while for galactic scales there 
is significant enhancement of the mass as detected from far away. 
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